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1.0  INTRODUCTION 


Ocean  Prediction  System  (OPS)  is  under  development  to  provide  estimates  of  ocean  surface  cur¬ 
rent  fields  and  surface  wind  fields  for  use  by  search  planners  in  real  time. 

These  estimates  will  be  used  by  CASP  2.0  to  develop  search  plans.  These  searches  should  be 
more  efficient  and  more  effective  because  in-situ  small  scale  dynamic  data  will  be  used  instead  of 
large  scale  historical  data. 

Search  effectiveness  is  expected  to  increase  due  to: 

•  a  greater  likelihood  that  the  search  region  will  contain  the  search  object,  and 

•  a  reduction  in  the  growth  of  the  search  area. 

Reference  [a]  provides  a  framework  for  conducting  a  cost/benefit  analysis  of  Ocean  Prediction 
System  effectiveness.  However,  Reference  [a]  assumes  the  subject  of  the  search  survives.  Surviv¬ 
ability  can  significantly  affect  the  cost/benefit  analysis,  due  to  the  high  cost  of  a  human  life. 

This  report  presents  a  methodology  for  treating  survivability  in  the  cost/benefit  analysis.  This 
approach  can  be  used  in  developing  Section  2.1,  Aspects  of  Survival,  in  Reference  [b]. 

The  report  is  organized  in  the  following  manner.  Factors  affecting  mortality  are  discussed  in  Sec¬ 
tion  2.0.  In  Section  3.0,  the  concept  of  conditional  probability,  which  is  used  to  combine  depen¬ 
dent  events,  is  described.  A  statistical  model  that  can  be  used  to  estimate  survivability  in  SAR 
cases  is  outlined  in  Section  4.0.  An  example  calculation  of  probability  of  survival  is  given  in  Sec¬ 
tion  5.0.  In  Section  6.0,  we  describe  how  to  combine  probability  of  detection  with  probability  of 
survival  in  order  to  estimate  the  probability  of  detecting  survivors.  A  report  summary  and  recom¬ 
mendations  are  given  in  Sections  7.0  and  8.0,  respectively. 

This  report  is  preliminary  in  that  a  number  of  assumptions  were  made  in  order  to  facilitate  sample 
calculations  using  the  proposed  methodology.  These  assumptions  should  be  checked  with  data 
where  possible.  In  addition,  a  simplified  model  for  cumulative  probability  of  detection  was  used. 


2.0  FACTORS  AFFECTING  MORTALITY 

The  probability  of  a  person  surviving  an  incident  at  sea  depends  upon  a  variety  of  factors  (Refer¬ 
ence  [c]). 

Factors  include: 

type  of  incident, 

person’s  physical  characteristics, 

environmental  conditions, 

use  of  flotation  devices  and  protective  clothing, 

shark  attack,  and 

food  and  water  reserves. 
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Incident  type  can  be  used  to  determine  how  likely  the  person  is  to  survive  the  initial  incident.  It 
can  also  give  a  perspective  on  how  likely  the  person  is  to  have  survival  gear  and  whether  or  not 
the  subject  was  injured  during  the  incident.  Examples  of  incident  types  include  a  ship  sinking  due 
to  weather,  a  ship  fire,  and  a  person  falling  overboard. 

Physical  characteristics  such  as  age,  weight,  and  fitness  can  affect  how  long  a  person  could  sur¬ 
vive  an  incident  at  sea.  In  addition,  training,  background,  and  personality  can  affect  survivability 
and  can  also  be  modeled.  However,  these  factors  will  not  be  addressed  in  Sections  4.0  and  5.0  in 
order  to  simplify  the  discussion. 

Environmental  conditions  can  influence  the  likelihood  of  the  person  surviving  the  incident,  as 
well  as  the  chance  of  survival  over  time.  Environmental  factors  include  temperature  (air  and 
water),  wind,  sea  state,  and  precipitation. 

Different  types  of  flotation  aids  can  affect  the  probability  of  survival.  For  example,  the  subject 
may  be  floating  unaided,  clinging  to  floating  debris,  or  in  a  life  raft.  Use  of  these  aids  may  result 
in  different  probabilities  of  survival.  These  devices  may  also  alter  the  probability  of  detection  by 
increasing  visibility.  Heavy  clothing  or  a  survival  suit  can  increase  the  chance  of  survival  against 
hypothermia. 

Food  and  water  reserves  can  affect  survivability  if  the  person  survives  other  immediate  effects 
(such  as  hypothermia  and  shark  attacks)  until  a  time  at  which  thirst  and  hunger  become  critical 
factors.  Reserves  are  a  combination  of  the  provisions  available  at  the  time  of  the  incident,  plus 
any  additional  provisions  that  may  be  obtained,  such  as  fish  and  rainwater.  Survival  equipment 
can  be  used  to  augment  and  extend  food  and  water  reserves,  and  to  favorably  impact  other  surviv¬ 
ability  factors. 


At  water  temperatures  above  68  °F ,  shark  attack  becomes  a  significant  risk  factor.  This  assump¬ 
tion  is  supported  by  the  graph  in  Figure  1-3  of  Volume  ff  in  Reference  [c]  which  is  based  on  U.S. 
Navy  data. 


3.0  MATHEMATICAL  PRELIMINARIES 

Some  of  the  mortality  factors  depend  on  one  another.  Therefore,  a  statistical  concept,  conditional 
probability  that  can  be  used  to  analyze  these  dependent  events. 

3.1  Conditional  Probability 

Conditional  probability  is  used  to  compute  the  probability  of  events  that  are  dependent  on  each 
other.  The  conditional  probability  of  event  A  occurring,  given  that  event  B  occurs  is  denoted  by 

P(A\B).  (3.1) 
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For  an  event  A  and  a  discrete  or  continuous  collection  of  possible  scenarios,  S ,  affecting  A ,  the 
probability  of  A  occurring  is  given  by  equations  (3.2)  and  (3.3),  respectively. 

For  a  discrete  set  of  scenarios. 


P(A)  =  5;P(A|5..)F(5,)  (3.2) 

i  =  1 

where  5,  are  the  different  scenarios  that  can  occur. 

For  a  continuous  set  of  scenarios. 


P{A) 


|p(A|S  =  s)dFs^,y 

R 


(3.3) 


=  Jp(A|5)JF5 


R 


(3.4) 


where  P{A\S  =  s)  is  the  probability  of  A  occurring  given  that  S  is  at  level  s  and  is 

the  probability  that  S  is  at  level  s  and  we  integrate  the  possible  range  R  of  values  of  s .  Equation 
(3.4)  is  a  commonly-used  shorthand  notation. 


3.2  An  Example 

Suppose  that  survival  only  depends  on  whether  the  person  is  wearing  a  life  jacket  or  not.  Let  L 
be  the  event  that  a  person  survives.  We  are  interested  in  the  probability,  P{L),  that  a  person  sur¬ 
vives.  We  condition  on  the  fact  that  a  person  is  either  wearing  a  life  jacket  J  or  not.  That  is,  the 
event  that  a  person  is  living  can  be  partitioned  into  the  event  that  a  person  is  living  and  wearing  a 
life  jacket  and  the  event  that  a  person  is  living  and  not  wearing  a  life  jacket.  Mathematically, 

L  =  {L  and  7}or{L  and  no  7}  .  (3.5) 

The  two  events  on  the  right  side  of  the  equality  are  mutually  exclusive  in  the  sense  that  one  cannot 
be  both  wearing  a  life  jacket  and  not  wearing  life  jacket  at  the  same  time.  This  means  that  the 
probability  of  the  left  hand  side  is  the  sum  of  the  probabilities  on  the  right  hand  side: 

P{L)  =  P{L  and  /}  +  P{L  and  no  J} .  (3.6) 
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Now  further  assume  that  we  know  that  a  person  has  a  20%  chance  of  wearing  a  life  jacket.  We 
also  assume  that  if  the  person  is  wearing  a  life  jacket,  then  he  or  she  will  survive  95%  of  the  time 
and  that  if  the  person  is  not  wearing  a  lifejacket,  then  he  or  she  will  survive  only  25%  of  the  time. 

We  can  combine  this  information  using  the  statistical  method  of  conditioning.  Conditioning  says 
that  the  likelihood  of  an  event  and  a  scenario  is  the  likelihood  of  the  scenario  times  the  likelihood 
of  the  event  given  that  scenario.  In  our  example, 

PiL  and  J)  =  PiL\J)P(J).  (3.7) 

This  means  that  the  probability  that  a  person  survives  and  is  wearing  a  lifejacket  is  the  probability 
that  the  person  is  wearing  a  lifejacket,  P(J)  ,  times  the  probability  that  the  person  survives  when 
he  is  wearing  the  life  jacket.  The  scenario  is  whether  the  person  is  wearing  a  life  jacket  and  the 
event  is  whether  the  person  survives.  We  similarly  obtain  for  the  other  term 

P(L  and  =  PiUDiPiJ^)) .  (3.8) 


where  “  means  complement  or  “not.” 
Therefore, 


PiD  =  P(L  and  7) +  />(L  and  70  (3-9) 

=  P(L|7)P(7)  +  P(L|70P(70  (3-10) 

=  (0.95)(0.20)  +  (0.25)(0.80)  (3.11) 

=  0.19  +  0.20  (3.12) 

=  0.39.  (3.13) 


This  example  is  for  a  discrete  distribution,  because  the  scenario  has  only  two  possibilities:  either 
a  person  is  wearing  a  life  jacket  or  not.  If  the  scenario  has  an  infinite  number  of  possibilities  (e.g., 
water  temperature  is  in  a  range  from  50°F  to  60°F),  then  a  continuous  distribution  would  be 
used. 
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4.0  STATISTICAL  METHODOLOGY 


In  this  section,  we  describe  a  statistical  methodology  that  can  be  used  to  model  the  likelihood  of 
survival  based  on  the  factors  described  in  Section  2.0. 

We  consider  the  following  factors:  water  temperature,  shark  attack,  use  of  a  survival  suit  and  food 
and  water  reserves.  We  do  not  explicitly  model  the  effect  of  a  person’s  physical  characteristics  on 
survivability,  or  address  incident  types.  The  incident  considered  in  this  report  is  one  in  which  the 
person  survives  the  initial  incident  and  enters  the  water  without  injuries  (i.e.,  person  in  the  water). 

We  use  the  following  notation: 

L{t)  :  the  person  is  alive  at  time  t , 

T  :  water  temperature, 

S  :  the  person  has  a  survival  suit, 

A(t)  :  a  shark  attack  occurs  by  time  t ,  and 
R{t)  :  resource  (food  and  water)  reserves  at  time  t . 

We  use  the  following  probabilistic  notation: 

:  the  cumulative  probability  distribution  (cdf)  of  random  variable  X , 

:  the  probability  density  function  (pdf)  of  random  variable  X , 

P(E)  :  the  probability  of  event  E ,  and 

:  is  the  complement  of  event  E ,  where 

PiE^)  =  1-P{E)  .  (3.14) 

The  quantity  of  interest  is  the  probability  that  a  person  is  alive  at  time  t,  or  equivalently, 
P{L{t)). 

4.1  Water  Temperature  T 

To  estimate  the  effect  of  temperature  on  survivability,  we  express  the  probability  of  survival  as  a 
function  of  the  temperature  distribution,  and  of  the  conditional  probability  of  survival  at  a  given 
temperature.  We  condition  on  water  temperature  T  first,  because  some  of  the  other  events 
depend  upon  the  water  temperature  T  (e.g.,  shark  attacks  rarely  occur  in  cold  water). 
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Treating  T  as  a  continuous  random  variable  we  have  that. 


P(L(0)  =  \pW)\T)dFj 


(4.1) 


where  the  integration  is  over  the  range  of  the  random  variable  of  T . 

We  need  to  know  the  distribution  of  T.  We  can  assume  a  distribution  for  T  or  we  can  estimate  it. 

4.1.1  Estimation 

In  order  to  estimate  the  distribution  of  T,  we  need  data  on  the  sea  temperature  for  SAR  cases. 
Using  this  data,  we  can  estimate  the  distribution  parametrically  or  non-parametrically. 

We  can  also  use  historical  data  to  obtain  information  for  a  given  season  and  region  and  then  esti¬ 
mate  the  distribution  from  this  data.  For  example.  Table  1  presents  mean  sea  surface  tempera¬ 
tures  for  selected  seasons  and  locations  off  of  the  coastal  United  States.  This  data  was  extracted 
from  Naval  Oceanographic  Office  databases. 


Table  1.  Average  Sea  Surface  Temperature  (°F) 


Season 

Location 

Jan  -  Mar 

Apr  - Jun 

Jul  -  Sep 

Oct  -  Dec 

Cape  Cod  Area 

45 

50 

65 

52 

Virginia  Area 

66 

70 

80 

70 

Florida  Atlantic  Coast 

72 

76 

83 

78 

Gulf  of  Mexico 

73 

79 

85 

79 

Gulf  of  Alaska 

40 

43 

52 

46 

Southern  California  Area 

60 

58 

63 

63 

Washington  Area 

46 

50 

60 

55 

4.1.2  Examples  of  Continuous  Distributions  for  T 

If  a  continuous  function  for  the  temperature  distribution  is  available,  then  a  continuous  formula¬ 
tion  of  the  conditional  probability  should  be  used.  Possible  assumptions  on  the  distribution  for 

water  temperature  T  are: 


•  r  is  distributed  uniformly  within  ±  a  degrees  (F  or  C)  tolerance  of  its  mean  value  T 
for  a  given  region  and  season, 

r-Unif[f-a,7'  +  a]  .  (4.2) 
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In  this  case,  the  pdf  of  T  is 


dFj  = 


|l/(2a),  iflr-fUa, 

[O  ,  otherwise. 


(4.3) 


•  r  is  normally  distributed  about  its  region  and  season  mean  T  with  a  standard  deviation 

s, 

T-'NiT,s)  (4.4) 


In  this  case,  the  pdf  of  T  is 


dFr  = 


ijlii 


1  (T-T) 

2 


(4.5) 


for  any  real  T. 
We  then  have 


P(L(0)  =  jpiLit)\T)dFj.  . 


(4.6) 


4.1.3  Effect  of  Water  Temperature  on  Survivability 

Figure  1  (taken  from  Reference  [c])  presents  survival  time  as  a  function  of  water  temperature  for  a 
person  in  water  without  an  anti-exposure  suit.  Data  are  shown  for  three  body  types;  fast  cooler, 
average  cooler  and  slow  cooler. 

If  we  had  information  about  the  physical  characteristics  of  the  person,  we  would  use  this  informa¬ 
tion  here  and  only  use  the  curve  for  that  person’s  body  type  and  a  distribution  for  dying  before  this 
time  (e.g.,  a  shifted  exponential  with  some  instantaneous  death  rate). 

For  each  body  type  shown  in  Figure  1,  we  use  a  non-linear  least  squares  method  with  an  exponen¬ 
tial  model 


fiT)  =  +  (4.7) 

to  estimate  the  survival  time  curve  at  each  temperature,  T,  using  the  following  procedure:  At 
equally  spaced  temperature  intervals,  determine  where  each  curve  intersects  a  given  temperature. 
This  creates  a  data  point  for  each  curve.  Make  a  data  set  for  each  curve  at  several  temperatures  (at 
least  10)  and  then  estimate  the  curves. 
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Figure  1.  Survival  Time  versus  Water  Temperature  (Figure  4-1  p.  4  -  5  Reference  [c]) 


Once  the  curves  are  estimated,  label  them  as  ,  and  for  fast,  average,  and  slow 

coolers,  respectively.  For  all  T,  the  curves  are  ordered  as  <fa{j)  *^fs(T)  ■ 

For  a  given  water  temperature,  T,  the  set  of  three  curves  provides  a  range  of  survival  times. 

From  Figure  1,  we  see  that: 

•  It  is  highly  likely  that  a  person  will  survive  up  to  time  , 

•  It  is  highly  unlikely  that  a  person  will  survive  longer  than  time  . 

•  The  probability  that  a  person  will  survive  up  to  time  t  in  between  the  two  bounds  decreases 
as  we  move  from  to  . 

Because  the  curve  for  the  average  cooler,  /^(jO  ,  is  not  perfectly  in  between  the  other  two  bound¬ 
ing  curves,  a  uniform  decrease  for  survival  time  between  and  is  not  warranted.  How¬ 
ever,  we  assume  that  a  uniform  decrease  occurs  for  survival  time  between  the  two  pairs  of  curves 


That  is,  we  assume  that  L{t)\T  is  distributed  as 

0.5(Unif(/^(7^,/,(^)  +  Unif(4(^,4(7^)) .  (4.8) 

which  gives  the  following  probabilities. 


P(L{t)\S‘=,T) 


1 

if  ^-ff(T)’ 

i-o.5/~-^^y  , 

Ja(T)-JfiT) 

if  ff(T)  -  ^  -fa(T)’ 

0.5  0.5/"^^/  , 

Ja(T)  ~JsiT) 

if  fa(T)  -  ^  -fs(T)’ 

0 

VI 

P 

(4.9) 


We  can  use  this  assumption  to  obtain  the  distribution  for  L{t)\S,  T  by  making  the  assumption  that 
when  wearing  a  survival  suit,  one  will  survive  k  times  longer  than  if  one  were  not  wearing  a  sur¬ 
vival  suit.  A  typical  range  for  k  is  between  2  and  10  (p.  4  -  6  Vol.  I  Reference  [c]).  This  means 
that 


L(0|5,  T~  0.5(Unif(k4^,  kf^^^))  +  Unif(k4(r),  kf^^j^))  (4.10) 
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which  gives  the  following  probabilities, 


P(L(0|5^r) 


1 

9 

if 

ff(T)  -  -faiT)’ 

1- 

0  5  ^ 

Ja{T)  -JfiT) 

0.5 

U.  J  r  X  9 

Ja{T)  -Js(T) 

0 

9 

if  fs{T)  ^ 

(4.11) 


4.2  Shark  Attack 

We  next  determine  the  conditional  probability  of  whether  or  not  there  is  a  shark  attack  by  time  t . 


That  is, 

P{L{t)\T)  =  P{L{t)\A{t),  r)P(A(r)|r)  +  P(L(r)|A"(f),  T)P{A<=\T)  (4.12) 

We  assume  that  if  there  is  a  shark  attack,  the  person  does  not  survive.  That  is, 

P(L(0|A(0,7’)  =  P(L(0|A(0)  =  0.  (4.13) 

Using  this  information, 

P{W)\T)  =  P(L(0|A(0,P)P(A(0|7’)  +  ^(i'(0|A^(0,  WA‘^(0|7’)  (4.14) 

=  P{L{t)\A<={t),T)P{A%t)\T)  (4.15) 

which  gives 

P(L(0)  =  \piL{t)\T)dFj  (4.16) 

=  ^P{L{t)\A^{t),T)P{A^{t)\T)dFj.  (4.17) 
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4.2.1  Distribution  for  A  (0 


Assumptions  are  needed  in  order  to  determine  P(A{t)\T)  .  Because  sharks  are  unlikely  to  be  in 
cold  water,  we  condition  on  water  temperature  T.  We  assume  an  exponential  distribution  on 
A(r);  that  is, 

P{A(t)\T)  =  l-e-'(p(7))  (4.18) 


where 

p(T)  is  the  instantaneous  probability  of  a  shark  attack  at  water  temperature  T. 

Use  of  an  exponential  distribution  means  that  sharks  are  equally  likely  to  appear  at  any  time  and 
that  one  can  expect  a  shark  attack  every  l/(j)(T))  time  periods. 

p(T)  can  be  continuous  or  discontinuous  (that  is,  having  a  jump  discontinuity  at  a  temperature 
T  where  it  is  known  that  sharks  are  not  present).  We  assume  that  p(T)  =  0  for  water  tempera¬ 
tures  below  some  temperature  7a  .  Note  that  7a  and  p{T)  may  change  depending  on  the  region 
and  season. 

The  distribution  of  A(t)  can  be  estimated  by  knowing  when  sharks  are  in  a  given  region  and  how 
likely  they  are  to  attack. 

4.3  Survival  Suit 

We  next  condition  on  whether  the  person  has  a  survival  suit  (i.e.,  event  S).  We  assume  that  hav¬ 
ing  a  survival  suit  is  independent  of  water  temperature  and  animal  attack;  that  is,  you  either  have 
one  or  you  don’t. 

We  also  assume  for  simplicity  that  the  probability  of  having  a  survival  suit  P{S)  is  constant.  We 
then  have 

PiL(t)\T,A‘^it))  =  P(L(0|5,  T,A<=it))PiS\T,A<^{t))  + 

P(L(0|5^  7’,A‘^(0)P(5^|7’,  A‘^(0)  (4.19) 
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This  gives 


P(L(0)  =  \p{L{t)\T)dFj 


(4.20) 


=  \p{L{t)\T,A<^it))P{A%t)\T)dFj 


(4.21) 


L(0|5,r,A"(f) 


L{,t)\S^,T,A%t)  ^(50 


P{A%t)\T)dFj 


(4.22) 


We  need  distributions  for 

L{t)\S,T,A^{t)  and 
L{t)\S^,T,A%t)  . 

This  type  of  estimation  is  somewhat  complicated  for  the  following  reasons: 

The  data  (hours  of  survival)  for  a  given  temperature  and  body  type  are  both  right  and  left  cen¬ 
sored.  That  is,  if  a  person  is  rescued,  then  you  only  know  that  the  person  lived  at  least  that  long; 
the  person  could  have  lived  longer.  This  is  called  right  censoring.  Also,  if  you  recover  a  dead 
body,  then  the  person  died  before  this  time.  This  is  called  left  censoring.  Data  should  be  analyzed 
using  statistic^  techniques  based  upon  censored  data. 

It  is  unlikely  that  enough  observations  will  be  collected  at  a  specific  temperature  to  get  a  good 
estimate  of  the  distribution.  The  data  should  be  smoothed  over  a  range  of  temperatures. 

It  is  also  unlikely  that  enough  observations  for  an  exact  body  type  (because  it  is  continuous)  will 
be  collected  to  get  a  good  estimate  of  the  distribution.  The  data  should  be  smoothed  over  a  range 
of  body  types  or  people  should  be  grouped  into  a  few  body  types.  Training,  background,  and  per¬ 
sonality  factors  could  be  treated  in  a  similar  fashion. 

4.4  Food  and  Water  Reserves 

In  this  section,  we  describe  a  model  that  can  be  used  to  estimate  the  effect  of  food  and  water 
reserves  on  the  probability  of  survival. 

Let  P(0  be  the  amount  of  the  resource,  either  food  or  water,  at  time  t .  We  condition  on  whether 
we  have  the  resource  or  not. 
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That  is, 


PiLit)\S,T,A<=it))  =  PiLit),R(t)>0\S,T,A%t))+  (4.23) 

P(Lit),Rit)  =  0|5,  r,A"(0) 

=  P{Lit)\S,  T,A<^(t),R(t)>0)PiRit)>0\S,  T,A‘^{t))  + 

Pmt)\S,  T,A%t),R(t)  =  O)P(RiT)  =  0|5,  T,A‘^{t))  (4.24) 

We  need  to  determine  the  distributions  of: 

(L(0|5,r,A‘^(0,/?(0>0),  (4.25) 

(L(0|5,  T,  A^{t),  R{t)  =  0) ,  and  (4.26) 

(/f(t)|5,r,A^(0);  (4.27) 

and  for  P(L{t)\S*^,  T,  A^(/)) ,  the  distribution  of: 

(L(^)|5^^,A"(0,/^(0>0),  (4.28) 

(L(0|5^  T,A^it),Rit)  =  0),and  (4.29) 

(/?(0|5Sr,A‘^(0).  (4.30) 


Recall  that  conditioning  on  A‘^(/)  means  that  there  has  been  no  shark  attack  up  to  time  t . 

4.4.1  Assumptions 

In  order  to  determine  L(i)  |5,  T,  A^^it),  R(t)  >  0 ,  we  assume  that  if  you  have  some  of  the  resource 
at  time  t,  R(t)>0,  then  the  person  is  alive.  That  is, 

Pmt)\S,T,A^it),R(t)>0)  =  1.  (4.31) 

Note  that,  even  with  a  survival  suit  and  no  shark  attacks,  the  person  may  die  from  exposure  before 
the  reserves  are  exhausted.  So  the  conditioning  is  done  on  temperature  first  as  discussed  in  Sec¬ 
tion  4.1.  Assumption  (4.31)  is  valid  at  temperatures  for  which  P(L(t)|r  =1). 

For  L(0 15,  T,  A^{t),  R{t)  =  0 ,  we  assume  that  the  person  will  automatically  live  h  hours  before 
there  is  any  effect  due  to  not  having  the  resource;  that  is, 

P(L(0|5,7’,A^(0,i?(0  =  0)  =  1  for  t<h. 
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Note  that  h  depends  upon  such  factors  as  the  person’s  physical  condition,  water  temperature,  and 
whether  the  person  is  wearing  a  survival  suit  or  not. 

For  times  t>h,we  assume  that  there  is  an  instantaneous  probability  of  death,  1{R)  ,  due  to  the 
absence  of  the  resources.  That  is,  the  distribution  is  a  shifted  exponential  distribution  where: 


P(L(r)|5,r,A‘^(0,/f(0  =  0) 


1  ,  if  r</i,  and 

g(.-(t-h))liR)^  if  t>h. 


(4.32) 


Note  that  this  probability  only  depends  on  t.  So  that  if  t<h,  then  there  is  no  impact,  otherwise, 
we  weight  the  probability  by  . 

We  also  assume  that  after  H  hours  without  food,  the  probability  that  the  person  is  still  living  is 
some  small  probability  p/f  =  Pj^iS,  T) .  That  is, 

P(L(r  +  if)|5,  T,A%t),R{t)  =  0)  =  (4.33) 

This  means  that  P{,L{H)\S,T,A^{t),R{H)  =  0)=pR  .  Using  these  conditions  for 
L(t)\S,  T,  A^{t),  R{t)  =  0  determines  the  instantaneous  death  probability  1{R)  . 

To  determine  Rit)\S,  T,  A‘^{t)  ,  we  assume  that  there  are  no  reserves  at  time  0  and  there  is  no 
possibility  of  obtaining  more  reserves  after  time  0. 

This  means  that  PiR(t)  =  0)  =  1  and  P(R(t)>0)  =  0  for  all  r. 

We  now  have,  for  t<h: 

P(L(0)  =  lp{L(t)\T)dFT  (4-34) 

=  ^PiL(t)\T,A<^(t))P(A^it)\T)dFT  (4.35) 

=  |[P(L(015,  r,A^(0)^(-^)  +  ^(^(0|-5^  T,A‘^it))PiS‘^)]PiA^(t)\T)dFj.  (4.36) 
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=  j[P(L(0|5,  T,A^(it),R(t)  =  0))PiS)  + 

P(L(t)\S^,  T,A‘^at),  Rit)  =  0))P(50]P(A"(0|T)^Fr 

(4.37) 

and,  for  t>h: 

P(L(0)  =  e-('-*>^(«)jp(L(0)^Fr 

(4.38) 

=  e-('-*)^wJp(L(0|r,A^(0)P(A^(0|70^i^7- 

(4.39) 

=  e-(»-/')'wJ[P(L(0|5,  r,  A‘^(0)P(5)  + 

P(Lit)\S^,  T,A<^(t))Pm]PiA‘^it)\T)dFj  (4.40) 

=  e-0-*)'wJ[P(L(0|5,  T,A‘^{{t),Rit)  =  0))P(5)  + 

P(L(0|SS  T,A<^at),  R(t)  =  0))P(Sn]P(A<=it)\T)dFj  (4.41) 

A  more  realistic  model  for  determining  R{t)  would  be  to  assume  that  R(t)  is  a  non-homoge- 
neous  Poisson  process.  This  type  of  model  includes  the  following  features: 

•  The  reserves  start  at  a  given  level,  Rq  , 

•  For  each  region,  season,  and  water  temperature,  additional  reserves  (food  and  water)  arrive 
at  a  Poisson  rate  ,  and 

•  For  each  region,  season,  and  water  temperature,  the  subject  consumes  the  reserves  at  a  rate, 

,  that  can  depend  upon  having  a  survival  suit  and  specific  body  type.  For  example,  if  it 
is  cold,  you  eat  more  to  stay  warm.  However,  if  you  have  a  survival  suit,  the  suit  keeps  you 
warm;  so,  you  do  not  eat  as  much  as  you  would,  if  you  did  not  have  a  survival  suit. 

The  equations  for  this  type  of  model  are  not  presented  in  this  report. 
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5.0  PROBABILITY  OF  SURVIVAL:  AN  EXAMPLE 


In  this  section,  we  illustrate  the  combined  effect  of  the  models  presented  in  Section  4.0  using  an 
example. 

Figure  2  shows  probability  of  survival  curves  versus  time  for  a  person  in  water  without  a  survival 
suit  for  three  average  water  temperatures,  50, 65  and  80°  F. 

The  following  assumptions  were  made  in  calculating  these  probability  curves; 

•  The  person  is  alive  at  time  t  =  0. 

•  Water  temperature.  Water  temperature  is  uniformly  distributed  with  a  ±10°F  toler¬ 
ance. 

•  .Shark  attack.  Probability  of  shark  attack  is  0.5  per  day,  if  water  temperature  is  above 
68°  F,  and  0  if  water  temperature  is  below  68°  F. 

We  assume  that  the  instantaneous  probability  of  shark  attack,  p(T)  ,  is  constant  after  a  specified  tem¬ 
perature  =  68“F.  If  the  probability  of  a  shark  attack  in  a  24  hour  period  is  0.5,  then 

P(A(24)|r)  =  (5-1) 

=  0.5  (5.2) 


Solving  for  P(T)  ,  we  have  that 

p(T)  =  (-ln(0.5))/24  (5-3) 

=  0.03.  (5.4) 

So  the  instantaneous  probability  of  shark  attack  is  0.03. 

•  Reserves.  Reserves  have  no  effect  on  survivability  for  up  to  12  hours,  and  the  person 
dies  after  48  hours  without  food  or  water.  That  is,  h  =  12  hours,  H  =  48  hours  and 
=  0.001. 
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Probability  of  Survival 


EXAMPLE 


Figure  2.  An  Example  of  ProbabUity  of  Survival  versus  Time  for  a  Person  in  Water  at 
Three  Temperatures 
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6.0  PROBABILITY  OF  SUCCESS 

We  now  calculate  the  probability  of  success  (see  Section  2.H,  Reference  [b])  with  the  constraint 
that  the  person  be  alive  at  the  time  of  detection. 

The  probability  of  success  at  time  t  is  given  by 

POSit)  =  I  POD{t')POA{t')dt'  (6.1) 

where 

POD{t)  =  Pr{Detection  at  time  r'|person  is  in  the  search  area}  ,  (6.2) 

and 

POA{t)  =  Pr{Personisinthesearchareaattime  r'}  .  (6.3) 

6.1  An  Example 

We  illustrate  how  to  calculate  the  probability  of  success  in  detecting  a  survivor  using  an  example. 

We  assume  that  the  instantaneous  probability  of  detection,  POD{t)  ,  is  independent  of  the  proba¬ 
bility  of  surviving,  P{L{t)) .  That  is,  the  probability  of  detection  does  not  depend  on  whether  the 
person  is  alive.  Then  the  instantaneous  probability  of  detecting  a  survivor,  POD^{t)  ,  is  given  by 

POD^it)  =  POD{t-tQ)P{L{t))  ,  t>tQ.  (6.4) 

where  tg  is  the  time  at  which  the  search  begins.  We  assume  for  simplicity  that  fg  =  ®  ’  that  is, 
the  search  begins  at  the  time  of  the  incident. 

The  probability  of  success  in  detecting  a  survivor  by  time  t  is  given  by 

POS^it)  =  j  POD{t')P(L{t'))POA{t')dt  (6.5) 

*^0 

We  now  assume,  for  simplicity,  that  POA{t)  =  1  ;  that  is,  the  person  is  in  the  search  area. 

We  use,  as  an  example,  a  search  conducted  in  an  area  that  contains  the  person. 
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If  the  incident  position  uncertainty  is  given  by  a  bivariate  normal  probability  density  function, 
then  the  probability  of  success  as  a  function  of  time  for  a  random  search  is  given  by 


POS{t)  =  1-e-z 


(6.6) 


where 


VWt 

^  “  2jta^a/ 


(6.7) 


V  is  the  search  speed,  and  W  is  the  sweep  width.  The  instantaneous  probability  of  detection  is 
given  by 


POD{t)  =  = 


VW 


(6.8) 


6.2  Sample  Calculation 

We  illustrate  calculation  of  the  probability  of  detecting  a  survivor  using  the  following  example. 

We  assume  that  a  person  is  in  the  water  without  a  survival  suit  and  that  the  person  enters  the  water 
at  the  time  that  the  search  begins.  For  simplicity,  we  ignore  the  effects  of  shark  attacks  and  food 
reserves. 

For  an  HH-65  search,  we  assume  that  search  speed,  V,  is  90  knots,  and  sweep  width,  W,  is 
0.1  nm.  Also,  we  assume  that  the  search  area  uncertainty  is  given  by  a  bivariate  normal  distribu¬ 
tion  with  =  Gy  =  lOnm,  and  that  we  are  searching  in  the  correct  area. 

Suppose  also  that  SAR  cases  are  evenly  distributed  in  the  January  to  March  time  period  across  the 
three  Atlantic  Coast  regions  shown  in  Table  1.  Thus,  P(T)  =  0.33  for  T  in  {45,66,72}.  The 
non-linear  least  squares  models  for  the  three  hypothermia  fatality  curves,  corresponding  to  fast, 
average  and  slow  coolers  as  shown  in  Figure  1,  are  given  by 


(6.9) 

(6.10) 


=  0.14  +  0.030e0078r^ 

fa(T)  =  0.75 +  0.026e0  089r^  and 

4(D  =  0.97+0.035e™r. 


(6.11) 


The  temperature  set  is  {45, 66, 72}.  Using  equations  (6.9),  (6.10),  and  (6.11),  we  have  the  follow¬ 
ing  values: 


•^(45) 

=  1.14, 

/a(45)  “ 

fsi45)  -  ’ 

•^(66) 

=  5.30, 

fai66)  -  ’ 

fsi66)  ~  25.1,  and 

•^(72) 

=  8.38, 

fain)  -  ’ 

/s(72)  =  • 

(6.12) 

Figure  3  shows  P(L(f)),  the  probability  of  survival  until  time  t,  using  the  expression  for 
P(L{t)  I T)  given  by  equation  (4.9). 

The  Probability  of  Success  Curves  for  POSit)  and  POS^(t)  are  shown  in  Figure  4  where 
POS(t)  is  given  by  equation  (6.6)  and  POS^it)  is  given  by  equation  (6.5)  with  POD  given  by 
equation  (6.8),  P(,L{t))  given  by  equations  (4.9)  and  (6.12)  and  POA  =  1 . 

Figure  4  illustrates  certain  effects  that  result  from  including  mortality: 

•  The  two  probability  curves  POS  and  POS^  are  the  same  for  t  <  5  hours  because  the  prob¬ 
ability  of  surviving  that  long  is  high, 

•  There  is  little  increase  in  probability  of  successfully  detecting  survivors  after  time  t  =  20 
hours  because  the  probability  of  surviving  is  very  small  by  that  time,  and 

•  The  POS  curve  continues  to  increase  because  it  doesn’t  include  the  effect  of  mortality. 

7.0  SUMMARY 

A  methodology  has  been  presented  that  can  be  used  to  estimate  the  probability  of  survival, 
P{L{t)). 

Estimates  for  P{L{t))  can  then  be  used  with  a  model  for  probability  of  success  to  estimate  the 
probability  of  detecting  survivors  in  SAR  cases  for  use  in  the  OPS  cost-benefit  analysis. 
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Figure  3.  An  Example  of  Probability  of  Survival  versus  Time 
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Probability  of  Success 
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8.0  RECOMMENDATIONS 


We  recommend  that: 

•  The  assumptions  made  in  this  report  in  developing  the  model  for  P{L{t))  should  be  veri¬ 
fied  with  data, 

•  Data  should  be  collected  in  SAR  cases  to  assist  in  estimating  these  parameters  (see  Appen¬ 
dix  A),  and 

•  The  POS  model  described  in  Section  6.0  should  be  extended  to  include  the  effects  of  the 
time  difference  between  the  incident  occurrence  and  the  start  of  the  search  and  the  move¬ 
ment  and/or  growth  in  search  area  due  to  location  uncertainty  and  motion.  Such  a  model 
should  be  used  to  estimate  OPS  impact  on  search  effectiveness. 
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[BLANK] 


APPENDIX  A 


The  following  SARMIS  data  should  be  collected  for  each  SAR  case  and  included  in  the  CASP 
database  to  assist  in  modeling  search  and  rescue  effecctiveness. 

A.1  Survival 

time  of  incident 

type  of  incident 

initial  location  of  incident 

location  at  time  of  rescue/recovery 

temperature  (air  and  water)  -  estimated  at  intervals 

wind  and  sea  conditions  -  estimated  at  intervals 

use  of  flotation  devices  and  survival  suits 

availability  of  survival  equipment 

time/cause  of  death 

time  of  rescue/recovery 

physical  condition  at  time  of  rescue 

presence  of  sharks 

subject’s  survival  training  and  environment  experience  level  (high,  medium,  low) 

A.2  Search 

search  assets  -  number  and  type  (include  speed,  altitude,  even  fatigue) 
time  late  -  time  between  incident  and  search  assets  on  scene 
time  of  detection 
time  of  rescue/recovery 

search  probability  area  size  and  location  as  function  of  time 
estimated  sweep  width 
actual  search  track 

visibility  in  search  areas  -  estimated  at  intervals 


A-l 


